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Abstract. The nature of some baryonic resonances is still an unresolved issue. The case of the A^*(1535) is 
particularly interesting in this respect due to the nearby r/N threshold and interference with the A'^*(1650). 
The A'^*(1535) has been described as a threshold effect, as a genuine 3-quark resonance, or as dynamically 
generated from the interaction of the octet of baryons with the octet of mesons. In the scheme of dynamical 
generation, predictions for the interaction of the TV* (1535) with the photon can be made. In this study, we 
simultaneously analyze the role of the A*'* (1535) in the tvN — > nN and 7 A tvN reactions and compare 
to the respective amplitudes from partial wave analyses. This test is very sensitive to the meson-baryon 
components of the A'^*(1535). 



PACS. 25.20.Lj Photoproduction reactions - 13.60.Le Meson production 
actions - 14.20.Gk Baryon resonances with S=0 
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1 Introduction 

The unitary extensions of chiral perturbation theory 
{U xPT) have brought a new light to the meson-baryon in- 
teraetion, showing that some well-known resonanees qual- 
ify as being dynamieally generated. In this pieture, the 
Bethe-Salpeter resummation of elementary interactions, 
derived from chiral Lagrangians, guarantees unitarity and 
leads, at the same time, to genuine non-perturbative phe- 
nomena such as poles of the scattering amplitude in the 
complex plane of the invariant scattering energy z, which 
can be identified with resonances. Coupled channel dy- 
namics plays an essential role in this scheme, with the 
chiral Lagrangians providing the corresponding interac- 
tions of the multiplets; even physically closed channels 
contribute as intermediate virtual states. 

The yl(1405) and the A^*(1535) have been explained as 
meson-baryon {MB) quasibound states [H 2 , 3 , 4 , 5 , 6][7l[8] 
from the interaction of the meson octet of the pion (M) 
with the baryon octet of the nucleon (B), mediated by the 
Weinberg- Tomozawa term. It is interesting to note that, 
even without chiral Lagrangians, the use of basic interac- 
tions for the coupled channels calls for an interpretation of 
some resonances like the yl(1405) as quasibound states of 
the scattering problem [^ fTUlITT] . In the last few years, the 
unitary schemes have been extended to the low lying 3/2~ 



baryonic resonances |12ill3j and other quantum numbers 

mm. 
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Dynamically generated states can also appear in other 
approaches. At higher energies, the dynamics contained 
in meson exchange is important. Meson exchange mod- 
els of the ttN interaction with the analyticity property, 
developed over the years [l6lfT7lfT8[|19[[20] . deliver a pre- 
cise description of the partial waves in ttN scattering. In 
Ref. [12], the Roper resonance could be described as dy- 
namically generated, i.e., without the need of a genuine 
(3-quark) resonance state and, in a speed plot, resonance 
parameters were extracted. In recent studies pTll^ . the 
amplitude could be analytically continued to the complex 
plane and, indeed, poles of the Roper were found. 

Electromagnetic properties provide additional infor- 
mation about the structure of strongly interacting sys- 
tems. They allow for independent tests of hadronic mod- 
els: In the picture of dynamical generation, the mesons 
and baryons can form a resonance. As the photocouplings 
to these components are well known, predictions can be 
made that can be compared to experimental results. In 
case of the A^*(1535), photoproduction of 77 on p and n 
has been tested in Ref. [53]. There, the dependence 
of the A^^'^ and S^^^ helicity amplitudes has been evalu- 
ated. While there is a qualitative agreement with the data, 
the faster theoretical fall-off with indicates a need for 
additional ingredients in the model such as a small gen- 
uine 3-quark state. Similar conclusions have been drawn 
recently in Ref. [24] . 
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Many tests of hadronic models with eleetromagnetic 
probes have been done at the level of cross sections only 
(see, e.g. Ref. [26]), mainly because of the scarcity of data. 
One can also test the radiative decay of the resonance [271 
[55] , but its experimental extraction is tied to large uncer- 
tainties. A first step towards a test of the amplitude itself 
and its phase has been done in Ref. [23], and the relative 
signs between A^^^, ^1/2 ^'^'^ "^i/s could be tested which 
are in agreement with the results extracted experimen- 
tally. Yet, the phase of the photoproduction amplitude is 
still not accessed through such a study. 

However, a direct access to the photoproduction ampli- 
tude is indeed possible via the corresponding multipolcs, 
in this case, £'o+- These quantities can be extracted in 
partial wave analyses that also take into account the po- 
larization observables [29ll30| . This allows to access the 
relative phases of the N*{1535) in ttN and ttN — > 

ttN. As we shall show, this provides a very sensitive test 
for hadronic models of dynamical generation. In practice, 
however, one has still to resort to phenomenological anal- 
yses which are, in principle, model-dependent because one 
cannot access the amplitude directly from the existing ob- 
servables, as there is not yet a "complete" experiment 
available. Thus, the aim of the present study is to analyze 
the reactions ttN ttN and 7iV ttN at the amplitude 
level and compare to the phcnomcnologically extracted 
multipoles [^[50] . 

Observable quantities that allow one to distinguish be- 
tween hadronic molecules and more elementary states are 
urgently called for. For S'-wave states close to thresholds 
this is discussed in Refs. [31]. For a model test, a quantita- 
tive description of the data is important. This was already 
stressed in Ref. [321133] . where tttt and ttN scattering was 
analyzed within the chiral unitary approach and the in- 
terplay of genuine and dynamically generated resonances 
was investigated. 

For a quantitative data description of the S'-wave pion 
production, the A^*(1650) has to be included. It has the 
same quantum numbers as the TV* (1535) and lies close-by; 
the interference among resonances with the same quan- 
tum numbers plays an important role as has been recently 
pointed out in Ref. [H]. The 7V*(1650) substantially mod- 
ifies the position and residue of the N*{1535). Such an 
interplay of resonances may even lead to the disappear- 
ance of a resonance as seen, e.g., in the D13 partial wave 
discussed in Ref. [22\ 

For the dynamical generation of the iV*(1535), we stu- 
dy only the S'-wave pion production as provided by the 
Weinberg- Tomozawa term. This effective interaction is clo- 
sely tied to p meson exchange which, in its full dynami- 
cal treatment, contributes to all partial waves and not 
only to the S'-wave. Thus, the strong couplings to the KA 
and KIJ channels, which appear in the description of the 
A^*(1535) as a dynamically generated resonances, will nec- 
essarily have an impact on higher partial waves. While this 
issue can be studied in the framework of meson exchange 
models [33], it is beyond the scope of the present S'-wave 
model. 



2 Formalism 

2.1 Meson-baryon interaction 

Various approaches have been followed in the past to de- 
scribe the properties of the A^*(1535) in terms of mesons 
and baryons rather than by the quark degrees of freedom; 
in Ref. [T112] the chiral meson-baryon Lagrangian through 
NLO provides the driving interaction. In Ref. [6], an SU(6) 
fine-splitting in the chiral Lagrangian provides the inter- 
action. We follow here the work of Ref. [5] which relies on 
the use of the leading order Weinberg- Tomozawa term. 

The strength of the meson-baryon interaction is fixed 
by the pion decay constant In S'-wave projection, the 
interaction to leading order is given by 

^fW--a.^(fc° + fc^) 

x,/M±I^/M±Z (1) 

V 2M, Y 2Mj ^ ' 

with the channel indices the baryon mass Af, meson 
energies meson decay constants fi (/,rj frj, Ik), and 
the center of mass energy z. Note the explicit SU(3) break- 
ing from the values of Jk = 1.22 /^r and = 1.3 /^r. The 
value of /k has been recently evaluated to a higher preci- 
sion with the sfightly changed result of fx — 1.193 [35l 
l36] . The coefficient is the coupfing strength of the me- 
son and baryon, which is determined by the SU(3) group 
structure of the channel. For strangeness, S = 0, and 
isospin,/ = 1/2, the values of the coefficient Cij can be 
found in Refs. [5ll37] for the two net charge states Q = 0, 1. 
We work here in the particle basis rather than in the 
isospin basis which allows to take into account threshold 
effects from different pion masses in 7r° photoproduction 
[cf. Fig. [8]. The six channels for Q = 0, 1, respectively, 
are given in the first column of Table [1] In Eq. ([T]) , terms 
of the order p/M have been neglected and we will con- 
struct the phototransition amplitudes in agreement with 
this approximation. 

The amplitude in Eq. ([1]) is the input for the Bethe- 
Salpeter equation 

rpNP _ yNP _^ yNP Q rpNP 

= (1-T/^PG)-VNP. (2) 

The notation l/^P ^q. ([J) and T^P in Eq. ^ is mo- 
tivated by the usual decomposition of the amplitude into 
the pole (TP) and non-pole (T^^) parts as defined by Eqs. 
([HIS]) in the following section. The potential V^^ is factor- 
ized on-shell as described in detail in Ref. [32]. Then, the 
propagator G factorizes and can be explicitly evaluated 
with the result given, e.g., in Eq. (23) of Ref. [13]. The 
factorized propagator G is logarithmically divergent and 
evaluated with dimensional regularization. For a thorough 
discussion on the connection of this ansatz to the N/D 
method [32] and the relation to the full, non-factorizcd 
Bethe-Salpeter equation, see also Ref. [27ll33]. 

In contrast to the original work on the A^*(1535) [5], 
in the present study the ttttN channel is not taken into 
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account. The main reason is that a consistent implemen- 
tation of photoproduction in this channel is beyond the 
scope of this work; second, the parameterization of the 
ttN ttttN transition potential from Ref. [5] is only valid 
at lower energies of the energy range considered in the 
present study. From a practical point of view, the ttttA^ 
channel has been found to be small in the isospin 1/2 
channel |S] and neglecting it seems to be a valid approxi- 
mation. We may notice the absence of the multipion states 
in the isospin 7 = 3/2 channel. Also, in ttN r/N, the 
ttttN channels is necessary to reduce the ttN rjN cross 
section from 3 mb to the experimental value of 2.5 mb as 
pointed out in Ref. [20] . 

It should be pointed out that there are four subtraction 
constants a^w, arjN, cika, and aKS coming from the log- 
arithmic divergence in the factorized propagator G. These 
constants are regarded as free parameters and used to fit 
the experimental data. The values of the subtraction con- 
stants can be related to the regularization scale to be 
fixed at "natural" values of around —2. In Ref. [8], approx- 
imate crossing symmetry of the amplitude T is imposed 
by the renormalization condition T{y^ ^ fi) ^ ^(m) 
and the renormalization scale is related to the baryon 
masses |25j . A similar "natural" renormalization condition 
has been found recently in Ref. [24] . 

As strict crossing symmetry is violated at other ener- 
gies than the renormalization point anyways, we allow for 
deviations from the natural values of the subtraction con- 
stants. Interestingly, it has been pointed out in Ref. [53] 
that those deviations can be related to genuine, three- 
quark like propagator terms in the potential. Thus, as 
in the present approach such deviations are needed phe- 
nomenologically, this is a strong hint that some genuine 
quark state is necessary for a quantitative description of 
the A^*(1535) properties in the reactions ttN and 

TriV ttN studied here. See also Ref. [23] where some 
need for such components has been found in the study of 
the helicity amplitudes within a similar framework as the 
present one. 

Second, the hadronic interaction in the present ap- 
proach is mediated by the Weinberg- Tomozawa term, and 
the pole diagrams (terms with D and F) that are present 
at the same lowest order in the chiral expansion [2, are 
neglected. In the present approach, we consider contribu- 
tions from those terms, and also from the contributions of 
higher order Lagrangians which are not considered here, 
to be absorbed in the values of the subtraction constants. 

In Ref. [5], the four subtraction constants are fitted 
to the Sii and 5*31 ttN partial wave amplitudes and a 
pole in 5*11 is found at around 1537 — 37 i MeV. It should 
be stressed that the TV* (1535) is dynamically generated 
even when the natural values for the subtraction constants 
of Ref. [24] are used. In this case, we find the pole at 
zq = 1653 — 145 i MeV with large couplings into KS and 
rjN, in qualitative agreement with the fine tuned model 
of Ref. [5|. Thus, the N*{1535) and its coupling pattern 
preexist in the "natural" scheme of Ref. [24|, while the 
fine tuning of the subtraction constants done in Ref. [S] 
merely moves the pole a 100 MeV down in energy into the 



N*{1535) region. Due to its large couplings to KA and 
KS, the TV* (1535) is sometimes referred to as a KA, KS 
quasibound state. 



2.2 Genuine resonance states 

Eq. ([2]) allows for the formation of poles in the T matrix 
because the matrix 1 — V^^G can become singular. These 
poles in the complex plane of the scattering energy z can 
be identified with resonances on the physical axis (Im z = 
0). In the present approach, this concept of dynamical 
generation of the TV* (1535) is tested by also allowing for a 
genuine TV* (1535) via a pole term of the form l/(z — Mr) 
in the potential. Thus, when performing fits to different 
data sets, a genuine resonance can replace the dynamical 
TV* (1535) from the original fit of Ref. [5|. This would be 
a sign that the solution of Ref. [5] is unnatural, and that 
a genuine resonance state in the potential leads to a more 
natural data description. The criterion to distinguish both 
scenarios is the quality of the fits to the data. 

A second important point is the presence of the TV* (1650) 
in the Sn partial wave state; it can interfere with the 
TV* (1535), even far below the TV* (1650) nominal energy of 
z =1.65 GeV [22]. Resonance interference is a necessary 
ingredient for a realistic description, and thus, we allow 
in the fits for a second genuine resonance. Both genuine 
resonances of the form l/(z — Mr) acquire their respec- 
tive widths through the rescattering via the meson-baryon 
loops. 

For the S'-wave couplings of the genuine resonance 
states to the meson-baryon channels ttTV, ?/TV, KA, and 
KS, we choose a derivative coupling. Such a choice is 
necessary to fulfill the minimal requirements of chiral sym- 
metry and, in particular, to reduce the influence of reso- 
nances at the ttTV threshold. For the derivative coupling, 
we take the same energy dependence ~ as given in the 
Weinberg- Tomozawa term from Eq. ([1]) itself. In particu- 
lar, we factorize the N*MB vertex on-shell in the same 
way as the Weinberg- Tomozawa term, i.e. k'^ is the on- 
shell energy and not a loop integration variable. 

For the SU(3) couplings of the resonances, we assume 
the resonance to be a superposition of (anti)symmctric 
octet state (8), 8, antidecuplct 10 and 27-plet. This set 
of states can be uniquely mapped onto the four physical 
states ttTV, r]N, KA, and KS. 

Altogether, the modification of the interaction kernel 
from Eq. ([1]) by the presence of a genuine resonance is 
given by the sum of the non-pole part V^^ from the 
Weinberg- Tomozawa term and a pole part from the 
genuine resonance, 

T/NP ^V-. ^ yNP , yP _ yNP ,f.f)f 

V *J ~ V ^ I] ^ ij I J ^ J I 

k° 

= y- [58 Cs(j) +58' Ca(i) +5ToC<i(«) +527 C27(«)] 
D-^ = Z - Mr, (3) 

with the bare resonance masses Mr and the bare vertices 
Fi that depend on the meson energy fc*^ in channel i. D 
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Table 1. SU(3) coefScients for the resonance couplings to the 
meson-baryon channels in the particle basis. Both cases of net 
charge Q = and Q — +1 are listed. All values carry a square 
root, i.e. —3/20 means —^3/20. 
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Fig. 1. Decomposition of the hadronic amplitude in pole and 
non-pole part according to Eq. ((5|. Diagrammatic represen- 
tation of bare (dressed) vertices F (F), self energy E, bare 
(dressed) baryon propagator D (D) and the pole part T^. 



stands for the bare propagator. In Eq. ([3]), a factor of 
I/Ztt has been inserted to have dimensionless couplings g. 
In Table [1] the SU(3) coefficients c are listed; note that 
the baryon-first coupUng scheme is used throughout this 
study. These coefficients can be obtained from Ref. [38] , 
but in constructing the isospin eigenstates, the minus signs 
from the phase conventions for 7r+ and also have to be 
taken into account. These phases are necessary to ensure 
that the genuine states are purely isospin 7 = 1/2 states. 
The values of Cs , Ca, Cd, C27 quoted in Table[T]include these 
additional signs. 

It is convenient to decompose the full hadronic ampli- 
tude into a non-pole part and a pole part accord- 
ing to (suppressing the channel indices) 

T = V + VGT 

= T^+T^^. (4) 

The non-pole part is defined as the sum of those diagrams 
that do not have a 1-particle cut. This is identical to T^^ 
from Eq. ([2|). A short calculation shows that 



r = r^GT = (1 + GT^^G)r, 

r= + r= r , (5) 

where F and F are the bare and dressed vertices, re- 
spectively; indicates the transposed of the vector F 
in channel space. The vertices F, self energy S, and the 
pole part T^, are displayed in Fig. [1] 



In the case of more than one genuine resonance present, 
Eq. ([5]) becomes 

Srr' = rjGFr' ^ {I + G F^^ G) T,- , 

= (i + t^pg) n = (1-1/npg)"' n , (6) 

where the indices r and r' label the resonances; the sum- 
mation is over the resonances. 

With the formalism of this section, the hadronic part of 
the model is defined. There are altogether 14 free param- 
eters: 4 subtraction constants for the meson-baryon loops 
and 4 coupling strengths per genuine resonance state with 
bare mass Mr- The model allows for two genuine states 
and also the formation of dynamically generated poles. 
In particular, one of the genuine states can replace the 
dynamically generated A^*(1535) which will be a test of 
which description of the A^*(1535) is the more appropri- 
ate one. 



2.3 Photon couplings to mesons and baryons 

In this section we describe the photon couplings to the 
meson-baryon components defined in the last section. The 
photon couplings to the genuine resonance states will be 
discussed in Sec. 12.4.21 They are gauge invariant by them- 
selves and can be discussed separately. 

In principle, gauge invariance is ensured by systemat- 
ically and consistently coupling the photon to all mesons, 
baryons, and vertices in the rescattering series provided 
by the Bethc-Salpeter equation. This has been recently 
discussed [39] and realized in Ref. ^40j within the context 
of UxPT ■ In the present case, apart from working in a 
non-relativistic framework, we have a factorized version 
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Fig. 2. Coupling of the photon to loops in 7A'^ —> N* . This set 
is gauge invariant in the fully relativistic treatment. Diagrams 
(d), (e), and (f) vanish in the heavy baryon limit. 

of the Bcthc-Salpctcr equation in Eq. and it is pos- 
sible to determine subclasses of diagrams whose sum is 
transversal. 

We summarize at this point a discussion from Ref. 
P5] to provide a clear overview of what has been done 
in the literature [40ll23l [2| concerning the transversality of 
the phototransition amplitudes within the U xPT ■ A fully 
consistent treatment has been realized in Ref. [l^. We 
start the discussion with an approximation to this general 
scheme, realized in Ref. |23j . where several electromag- 
netic properties of the iV*(1535) have been evaluated. 

(I) : In Fig. [5] we show a set of diagrams with pho- 
ton couplings in the process — > N* with the quan- 
tum numbers of the TV* being 1/2^. For the discussion of 
transversality, we use the relativistic N and TV* propaga- 
tors and vertices. It has been shown in Ref. [23j that the 
setElJa) to[2][f) is gauge invariant. 

(II) : One can factorize the (hadronic) N*MB vertex in 
Fig.[2]on-shell, i.e., take the N* MB vertex out of the loop 
integration. In this case, it has been shown in Ref. [23] that 
diagrams (a) to (e) form a transversal set of diagrams. 

(III) : A further simplification of the hadronic part of 
the model is the heavy baryon approach as discussed in 
Ref. US] . In that case, the Dirac structure of the propaga- 
tors and vertices is simplified. Also in this non-relativistic 
formalism, the gauge invariance holds order by order in 
the 1/M expansion as discussed in Ref. [53]. Again, dia- 
grams dja) todl^e) form a transversal set of diagrams. In 
the non-relativistic treatment, one considers the convec- 
tion part of the -f{k)B{p) B{p') coupling, ~ e{p + p'), 
separately from the magnetic part, ^ k x p; the latter is 
transverse by itself. 

(IV) : A further simplification has been done, e.g., in 
Ref. |2]. In that case, even the phototransition amplitude 
itself is factorized on-shell, in addition to the on-shell fac- 
torization of the meson-baryon vertices: rather than eval- 
uating the full loops in Fig. [2l an effective range expansion 
of the tree-level diagrams in Fig. [3] is performed. The re- 
sulting phototransition amplitude is then multiplied with 
the factorized meson-baryon loop. We will investigate the 
consequences of such a procedure in Sec. 13.1.11 

The items (I)-(IV) above summarize what has been 
done in the literature. In the present study, we evaluate the 
phototransition amplitude at the level of approximation 
(III). This means we consider the photon loops (a) to (e) 



from Fig. [2] in the same non-relativistic framework that 
has been used for the hadronic part of the model discussed 
in Sec. 12. II In particular, we factorize the on-shell hadronic 
vertex of the type MB N* on the right-hand sides of 
the diagrams in Fig. [21 The rescattering scheme from Eq. 
([2|) provides also the Weinberg- Tomozawa term of the type 
MB ^ MB (cf. Fig.[l|). This vertex is treated on-sheU as 
well, i.e. it does not contribute in the loop integrations of 
the loops in Fig. [2] 

This scheme provides, in principle, a well-defined set 
of diagrams that can be evaluated straightforwardly. It 
can be realized in a field theoretical approach within an 
approximation as described in the next section. 

2.4 Formalism of photoproduction 

In the present investigation, we follow the work of Ref. [41] 
to construct a photoproduction amplitude based on the 
field-theoretical approach developed by Haberzettl |42| . 
The full formalism of Ref. [42 results in a complex and 
highly non-linear amplitude which requires truncation of 
some of the reaction mechanisms and/or replacement by 
phenomenological approximations for its practical imple- 
mentation. Any such an approximate treatment should 
preserve the relevant symmetries of the original full ampli- 
tude. In particular, unitarity and gauge invariance should 
be maintained. This has been done recently in Ref. |41j . 
where the complicated part of the interaction current has 
been approximated by a generalized contact current such 
as to preserve the gauge invariance (and unitarity) of the 
original full amplitude. 

In the following, we show how the treatment of Ref. 
[5T] is applied here. We denote the four-momenta of the 
initial state photon and nucleon by k and p, respectively. 
The four-momenta of the meson and baryon in the final 
state arc denoted by q and p' , respectively, s = {p + k)^ , 
u= (p' — fc)^ and t={q — ky. 

In this section we do not consider the genuine reso- 
nances; they will be included in Sec. 12.4.21 

The approximate photoproduction amplitude derived 
in Ref. [H] is expressed as 

M'^ = M^ + + M^ + A/f„, , (7) 

where M^ denotes the x-channel (x = s,u,t) tree-level 
amplitude (all involving physical masses and couplings) 
and Mj^^f , the interaction current. The latter is given by 

Knt ^ ^4" + T^ + T^^G [M^T + ^hr + T"] , (8) 

where M^^, (x ~ w, t) stands for the transverse part of 
M^, i.e., fc^M^y = 0. denotes the transverse contact 
current unconstrained by the Ward-Takahashi identity. As 
has been mentioned in [41] , it can be fixed from the data. 

(given below) denotes the generalized contact current 
which is an approximation to the complicated part of the 
interaction current that is not taken into account explic- 
itly 

The terms with G in Eq. ([5]) imply an integration over 
the loop four-momentum, which is not indicated explicitly. 
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(c) (d) 



Fig. 3. Tree-level contributions in photoproduction. Meson 
pole term (a), baryon pole term (b), KroU-Ruderman term (c), 
and crossed nucleon pole term (d). They correspond, respec- 
tively, to the t-channel Mj^ , s-channel Mif", KroU-Ruderman 
contact Af^^, and u-channel currents. 




Fig. 4. Photoproduction loop amplitude [second term in the 
r.h.s. of Eq. ([10} or Eq. J^l] constructed from M^j,, the 
transverse part of the t-channel diagram Fig. [3] (a), and from 
M^j,,the transverse part of the u-channel diagram Fig. [3](d). 
The intermediate loops imply a sum over allowed SU(3) states. 
The shaded circle indicates the factorized hadronic amplitude, 
T^-P from Eq. 

In particular, the loop function is not factorizing. In order 
to distinguish these terms from the previous formalism 
of the hadronic amplitude, where G is indeed factorizing, 
we denote the non-factorizing terms with G instead of 
G. Furthermore, throughout this paper, we refer to this 
loop (G) as the photon loop since it involves a photon 
attaching to one of the particles (meson or baryon) in 
the loop integral as illustrated diagrammatically in Fig. IH 
This is not to be confused with the intermediate jB loop 
which would be present in a coupled channel formalism 
beyond one-photon approximation. The present formalism 
is within an one-photon approximation. 

Unlike some meson exchange models, the present model 
contains no (phenomenological) form factors. In this case, 
the generalized contact current M^, as prescribed in Ref. 
[IT] , reduces to the usual (dressed) Kroll-Ruderman con- 
tact current 

M!^ = M^!^ ^ r^BBM, : (9) 

where r^gj^^^ denotes the BBMj Kroll-Ruderman con- 
tact vertex. 

As mentioned before, the transverse contact current 
T^, unconstrained by gauge invariance, can, in principle, 
be fixed from the data. In the present work, however, we 
set it to T'' = for simplicity; the existing data show no 
clear sign for its necessity. (See, however, the discussion of 
the £^0-1- multipole amplitude for neutral pion production 
in Sec. [321) 



The approximate amplitude [cf. Eqs. (|7|8p ] then can 
be written as 

Aff = V^+ T^^G [M^T + M/'y] , (10) 

where 

T/^ = + il/^ + il/f + Afj^^ (11) 

constitutes the sum of the usual tree-level Feynman dia- 
grams which is gauge invariant. It is illustrated diagram- 
matically in Fig.[3| The second term in the r.h.s. of Eq. pOj) 
is illustrated in Fig. [H 

The above equation defines the photoproduction am- 
plitude of the present model excluding the genuine reso- 
nances. Note that, in general, a regularization procedure 
is required in order to carry out the integral in Eq. (jTUj). 
In the remaining of this section, we show how this is done 
here. Using the notation 

M^{u{p')\M'^eMp)) . (12) 

where \u{p)) stands for the nucleon spinor without the 
Pauli spinor and the normalization {u{p)\u{p)) ~ 1, Eq. ([T0|) 
becomes, 

M^V + A~Iioop. (13) 

where 

Mioop = {u{p')\T''''G [Mi^j, + M/^j,] 6^|u(p)) . (14) 

In the present work, we calculate the integral in the above 
equation in the non-relativistic limit (|p| <C M). Also, our 
model uses an on-shell factorized final state interaction 
yNP \Yjtj^ these approximations, it can be shown [I51I231 
[55] that the loop amplitude Mioop in the above equation 
decomposes into the form (in the cm. frame of the system) 

Mioop = e^r^V, , (15) 

where cr,y = (0, cr) and 

r^'' = a g^" + b Pt^P" + c k^'P" + dPfk'' + ekfk'' , (16) 

with P = p + k = q + p' = [P" = z, 0). 

By construction, I^^" is transverse, i.e., k^F'^'^a^ = 
since Mioop is transverse. This leads to the relation among 
the coefficients in Eq. 

a + dk-P + ek^ = . (17) 

Note that in Eq. ([T6]) the only terms that actually con- 
tribute to the loop amplitude in Eq. ([T5]) arc the a and 
d terms. The other terms contribute nothing once they 
are contracted with ai, and the photon polarization vec- 
tor e^. Therefore, within the present approximation, the 
loop amplitude has the structure 

Mioop = —CL cr ■ e — d (7 ■ k Ef) P'^ . (18) 

For photoproduction, only the a term contributes in the 
above equation, while both the a and d terms contribute to 
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electroproduction. In the former reaction, from Eq. (|17p . 
the coefficient a is related to coefficient d by 



(19) 



-dk- P 



It should be noted that although the e term in Eq. (|16p 
does not contribute to the loop amplitude Mioop [cf. Eq. (fT8 
it is crucial for the proper gauge invariance condition given 
by Eq. (jlTp which is relevant for electroproduction. 

It is instructive to compare the structure of the loop 
amplitude given by Eq. (fT5|) with the most general form 
of the electroproduction amplitude [44] , 



iMfuu = iJicr ■ e + J2(t ■ qcr 
+ iJ^cr ■ qq ■ e — ijjcr 



(k X e 



iJscr ■ kq ■ 
iJscr ■ keQ . 



(20) 



A comparison of Eq. (jlSp with the above equation shows 
that the a term corresponds to the Ji term and the d term 
to Jg. All other structures J2 to J7 are zero in the loop 
amplitude Eq. It is also immediate that the .S-wave 
state can only contribute to the Ji and Js terms in the 
above equation (a and d terms in Eq. ([T8|)). 

As mentioned before, our model treats the loop inte- 
gral in Eq. (|14p in the non-relativistic limit, in addition 
to using an on-shell factorized T^^ . With these approx- 
imations, the coefficient d in Eq. p8)) can be extracted 
straightforwardly by a direct integration of the loop in 
Eq. (fH)) rendering a well defined (finite) value, a feature 
that can be shown from dimensional considerations. On 
the other hand, the direct loop integration yields infinity 
for the coefficient a which, as mentioned before, calls for a 
regularization of the loop integral. We do this effectively 
by determining a via Eq. (fT9|) . In the following subsection, 
some details on the evaluation of the coefficient d is given. 

Eq. pO|) can be extended to a coupled-channel ap- 
proach as 



Mt = vr 



uT 



(21) 



where the subscript i(j) specifies the baryon- meson chan- 
nel. 

Finally, before leaving this section, we mention that 
the present photoproduction amplitude given by Eqs. (jlOl 
lll[) can also be obtained from the full amplitude given by 
Eq. (3.20) of Ref. g^, provided we identify the term 
of Ref. IH] — which is to be approximated by a gener- 
alized contact current — with the sum of the bare KroU- 
Ruderman term (the first diagram in Fig. 3(F) ) and dia- 
gram Fig. 3(G) of Ref. [ID]. In addition, the dressed pho- 
ton vertex of Ref. [41] should be identified with the dressed 
photon vertex of Ref. 00], represented there in Fig. 3 by 
the photon line attached to the open square. 



2.4.1 Evaluation of the photon loops 

In this section we give some details on the calculation of 
the coefficient d appearing in Eqs. P^[TB|) . In the follow- 
ing, since we are interested only in the jN Ntt process 



Table 2. Factors 
charge channels j. 



Aj for the photon loops in the different 



For the channels K^E". K°A, 



and 



Vn, Aj 



■ 0. 





TT+n 




VP 






-(D + F) 


\/2(D4 


-F) 


-{3F- 


D)/V3 


D-F 


K+A 






K+E- 




■K~p 


-{D + 3F)/VS 


-V2{D 


-F) 


V2{D 


-F) 


-V2{D + F) 



in the present work, we restrict ourself to this channel, 
i.e., i = Ntt in Eq. ([2T|) . 

With the on-sheU factorization of the FSI, T^^, the 
amplitude from Eq. ([TBI) can be written as (i ~ Nn) 



4 1! 



d^p 



3 ' 



Mut{p") + Mtrip") G,{p") 
i 



(22) 



where we have displayed only the relevant argument on 
which each quantity depends upon. Fj can be expressed 
in the form 



Fj = -—Ajk-PdjCT-e, 



(23) 



where factors Aj for channel j combine the SU(3) factors 
of the MBB and the jMM vertices. They are given in 
Table[2]for each channel, dj is the part of the loop integral 
contributing to the coefficient d in Eq. ([TH]) expressed in 
terms of the Feynman parameter integrals. The contribu- 
tion to it from diagram |4K a) is finite with the result 



Ma) 



(47r)2 



1 1 
dx 



dz 







x{z - 1) 



x[{x - l)s + z{s - Mfj) + M2] + {1- x)m 



2 ' 



(24) 



where Mj (rrij) is the baryon (meson) mass in the loop, 
Afjv is the mass of the incoming nucleon (proton or neu- 
tron). The factor \ej\ is the modulus of the charge of the 
loop meson (0 or +e). 

The contribution to dj from diagram[4Ub) is subleading 
in 1/M. Yet, it is evaluated for a discussion on theoretical 
uncertainties in Sec. 13.31 with the result 



#^) = |e,,.|l^ /d. / dz 



(47r)2 





X z 



' x[{x - l)s + z{s - Ml) + m2] + (1 - x)M'j 



(25) 
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where \ eBj \ is the modulus of the charge of the loop baryon. 
We have used here only the convection part of the jBB 
coupling as discussed in the following. 

Diagram[4ljb) contributes to the d term but the contri- 
bution is subleading [5711^ in 1/M. We neglect this term 
for the numerical results, but consider it in the discussion 
on theoretical uncertainties in Sec. 13.31 

There are various contributions to diagrammb) as dis- 
cussed in detail in Ref. [23]. Apart from the convection 
part of the jBB coupling, evaluated in Eq. ((25|) . there is 
a magnetic part. Both parts are of similar and of small 
sizes [23]. Furthermore, in Ref. [23] the S'^A transition 
magnetic moment has been considered, whose contribu- 
tion also results to be small. In particular, all these differ- 
ent contributions for diagram IH^b) are not only small in 
size, but barely change the phase of the phototransition 
amplitude (cf. Table V of Ref. [23]). 



2.4.2 Photocouplings to genuine resonance states 

Apart from the phototransitions via loops, we also allow 
for direct jNN* coupling to the genuine resonance states. 
For the photon couplings to a genuine N* = Sn resonance 
we consider the Lagrangian 



C 



iV*75 



2M, 



(7u. d-" A'' N 



N 



(26) 



which provides a pure transversal 7iV N* transition. 
The transversality holds order by order in momentum, 
which can be also shown explicitly. In the non-relativistic 
reduction that is used in this study, Eq. pS]) leads to the 
vertex 



F - 



2Mn 



M 



(27) 



where k is the momentum of the incoming photon and 
gjNN* is a free parameter. In the present model, we allow 
for two genuine resonances, which both appear in charge 
zero and charge -1-1. Thus, the vertex from Eq. ([27]) ap- 
pears four times in the model, with four independent cou- 

pHng constants g^^p'^,, 4n^*- 

The inclusion of genuine resonance states leads to an 
additional contribution to the photoproduction amplitude 
given by Eq. ([7|. This additional contribution is given by 



1 



(28) 



where the summation runs over the resonances labelled by 
indices r and r'. In the present study, the N*N^ electro- 
magnetic transition vertex, Fj^,, in the above equation is 
dressed according to 



(29) 



for a given resonance labelled r'. A diagrammatic repre- 
sentation of Eqs. ([^I29p is shown in Fig. [5] In the actual 
calculation, the jNN* vertex from Eq. ([77]) is factorized 




Fig. 5. Photoproduction amplitude with s-channel genuine 
resonances, corresponding to Eq. (|28|l . F'^ (F'^) is the bare 
(dressed) photon vertex [Eq. (|29p ]. F (S) is the dressed 
hadronic vertex (self energy). 



on-shell. Note that, in principle, the genuine resonances 
can also contribute to the u-channel amplitude Af^. We 
have ignored this contribution in the present work. 

With the inclusion of the genuine resonances as de- 
scribed above, the full photoproduction amplitude in the 
present work becomes 



M 



full 



N* 



(30) 



where M is given by Eq. ([T3]) and 
2.4.3 Tree-level amplitudes 

Since in the present work the hadronic processes are de- 
scribed in the S partial- wave state [cf. Sees. 12.11 12. 2j . 
the nucleon s-channel amplitudes A4l^ won't be dressed 
by the hadronic interaction in the P-wave state. Further- 
more, since the loop integrals are evaluated in the non- 
relativistic limit, there are no dressings of the nucleon am- 
plitude in the S'-wave (negative energy) state either. 
Therefore, we evaluate the nucleon contribution to the 
tree-level amplitude V in Eq. ([13]) [cf. diagrams in Fig. [3])], 
including the s-channel amplitude A/^, using the physi- 
cal nucleon mass and physical coupling constants. This is 
also in accordance with the gauge invariant approach of 
Ref. [40] , where the photon-baryon vertex (represented by 
the photon line attached to open square in Fig. 3 in that 
reference) is not dressed by the equation. In contrast, the 
genuine 5ii resonances are dressed as mentioned in the 
previous section. 

Also, the nucleon tree-level amplitudes are evaluated 
with the full Dirac structure and relativistic propagators 
and then projected onto 5-wave to filter the S-wave con- 
tribution to photoproduction wc are interested in. 
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The photoproduction amplitude is constructed from Similarly, for the j{k)N MB process, the total cross 
i;ian densities 

TV (75 ^tvt)N , (31a) 



the following Lagrangian densities 
D + F 



sections are given by 

1 MnMb Qmb 



NN-n 



(J. 



2/. 



2M 



N 







D + F - 
e——N-i^4[T X 7r]3iV 



A^] N , (31b) 

(31c) 
(31d) 



Note that in the above Lagrangian we have also taken 
into account the anomalous term, which is higher order 
in photon momentum. Here, A'' and tt denote the nucleon 
and pion fields, respectively. The vector notation refers 
to the isospin space. We take D + F = 1.37. denotes 
the electromagnetic field, e is the proton charge; e = (1 + 
r^)/2 and k = [1.79(1 + r^)/2 - 1.93(1 - r^)/2] are the 
charge and anomalous magnetic moment operators of the 
nucleon, respectively. 

The propagators required for constructing the tree- 
level amplitudes are 



1 



S{P) = 



1 



i)-M, 



N 



(32) 



where A{q) denotes the pion propagator with mass 
and, S(p), the nucleon propagator with mass Mn- 



2.5 Observables and analysis 

Summarizing the present model, we have four free pa- 
rameters from the subtraction constants of the channels 
ttTV, TyiV, KA, KS (6 channels in the particle basis). For 
the two genuine resonance states allowed in the model we 
have, for each of them, four couplings to meson-baryon 
plus two couplings to the photon as free parameters, in 
addition to the bare mass. Altogether, we thus have 18 
free parameters. 

In order to compare to cross sections and/or partial 
wave analyses, we quote their connections with T and 
Mfuu as defined in this study. For ttN — > AIB, the hadronic 
amplitude T given in Sec. 12.21 is connected to the dimen- 
sionless partial wave amplitude T from Ref. [45| by 



TttN^MB 

Paib 



— ^/PttN Pmb Tt^N- 

Mb Qmb 
47r z 



>MB, 



(33) 



where Mb is the mass of the outgoing baryon and Qmb is 
the magnitude of the cm. three-momentum in the outgo- 
ing MB channel. The total cross sections for ttN MB 
are 

1 Qmb Mn Mb . ,2 

(^itN~*MB — -; — \-'-itN^MB\ ■ yo'i) 



m'mA 



|fe| 

[^m'\M^j^_MB\^m) 



, (35) 



where M^^ 

^MB stands for the photoproduction ampli- 
tude, Mfuu, given by Eq. ([30]) . A stands for the two inde- 
pendent polarizations of the photon beam, m! (m) denotes 
the spin-projection quantum number of the nucleon in the 
final(initial) state. 

Since we restrict ourselves to the S partial-wave only, 
the photoproduction amplitude Mfuii is related to the fa- 
miliar multipole amplitude i5o+ |44j by 



^Mn Mb , 

^0+ — — -. Ji- 

47r z 



(36) 



where Ji is from Eq. (|20p . 

In terms of £'0+ j the photoproduction total cross sec- 
tions are 



a. 



'-yN^BM 



Ifcl 



+ \jN^MB 



(37) 



The Eq+ multipole in the isospin basis for the outgoing 
particles is given by 

72 1 

pEo+{Sn) = ^ Eo+imT+) + - Eo+{pn°), 

\/2 1 

„£'o+(S'ii) = — Eq+{ptt^) - - Eo+{mT°), 

Eo+iSs,) = Eo+ipir'^) - 4^ Eo+ip7r+), (38) 

v2 

that show an additional factor of 1/ \/3 compared to Clebsch- 
Gordan coefficients. 

The T and Mfuu amplitudes have poles in the com- 
plex plane of the scattering energy z. These poles lie on 
the unphysical sheets. It is thus, necessary to analytically 
continue the amplitude to these sheets. This is a standard 
procedure and described in detail in Sec. 14.11 Poles of T or 
M full in the complex plane can come from the genuine res- 
onance states in the model, but also from the unitarization 
via the Bethe-Salpeter equation ([2]); indeed, the non-pole 
part from Eq. (j4]) itself obeys the Bethe-Salpeter equation. 
In particular, the term 1 — V^^G'-^'' can become singular 
if the Weinberg- Tomozawa term provides sufficient 
attraction. This will lead to a "dynamically generated" 
pole in T. 

To analyze the poles, it is convenient to perform a 
Laurent expansion at the pole position at complex zq- The 
leading term provides the residue a_i and the hadronic 
amplitude can be written as 



a-i 



a~i 

Z- Zq 

9t 9] 



+ ao + ai{z — zq) + 0(2:^), 



(39) 
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Fig. 6. Multipoles Eo+ in Sn. The red band shows the pre- 
diction using the full and simplified model from Ref. [5]. The 
solid red line shows the prediction of a refit to ttN scattering 
within the full model of Ref. [5]- The theoretical curves have 
been shifted upwards/downwards to have a better overlay with 
the data. Black data points: SES yN — » nN partial wave anal- 
ysis from Ref. [FA07]. Gray points: SES from Ref. [2S] 
[FA06]. Solid (dashed) black lines: Analyses SAID from [45] 
[FA07] (MAID2007 from [30]). 

for a transition from channel i to j. The residues are pa- 
rameterized as products of values gi gj which we call cou- 
pling strengths. We call Eq. ([39]) . through the residue term 
with a_i, the pole approximation (PA) of the amphtude. 



3 Results 

3.1 The i?o+ multipole and the phase problem 

In this section, we will see that there is a serious phase 
problem in Eq^ tied to the model of Ref. [S] . We will also 
see that in order to resolve the issue, we have to allow for 
the presence of additional resonances. 

Before coming to the results of the present model in- 
troduced in Sec. 12.11 we show the i?o+ multipole evalu- 
ated from the original model of Ref. [5] . In that model, no 
genuine resonances are present and the A^*(1535) is fully 
dynamically generated from the coupled channel interac- 
tion and the unitarization from Eq. This means, the 
interaction kernel V is entirely given by the Weinberg- 
Tomozawa term V^^ from Eq. ([1]). As there are no bare 
jNN* couplings, the photon interaction is given entirely 
by M from Eq. ([1^ . 

The results for pEo+ (photoproduction on the pro- 
ton) and nEo+ (neutron), using the model from Ref. [S], 
are shown in Fig. [H] To compare the energy dependence, 
we have slightly shifted the theoretical results upwards/ 
downwards to match them to the experimental i?o+ (in 
the following, we ignore this discrepancy and concentrate 
on the resonance shape) . The red band in Fig. [S] is given 



by two different solutions (full and simplified) from Ref. 

in the simplified version of the original full model the 
ttttN channel and the form factors from the Weinberg- 
Tomozawa term are omitted [5]. Yet, as Fig. [5] shows, the 
outcome is not very sensitive to these details. 

In the model of Ref. [5] , there are residual discrepan- 
cies with the ttN phase shifts. The deviations observed for 
i?o+ in Fig.[S]may be due to these discrepancies. Thus, we 
have performed a refit to ttN ttN using the full model 
from Ref. [5] , but choosing only a narrow energy interval 
around the A^*(1535) and not fitting to the 6*31 partial 
wave amplitude. The fit parameters are the four subtrac- 
tion constants for the loop functions of the channels 
ttN, 7]N, KA, and KS. In the refit the theoretical solu- 
tion in ttN — *■ ttN matches much better the results of the 
TriV PWA analyses [iM^ in the iV*(1535) region. In par- 
ticular, in the refit the A^*(1535) becomes wider while the 
fitted subtraction constants are still close to their original 
values quoted in Ref. [S]. Yet, when evaluating £^0+ from 
this refit, the result stays qualitatively the same, indicated 
with the red solid line in Fig. [6] 

Thus, the A/'* (1535) from the model of Ref. [5] is seri- 
ously off phase in i^o+j and this result does not depend on 
the details of that model, nor can it be easily cured with 
a refit, which delivers a better agreement in TriV ttN. 
Note that the Watson's theorem relates the phase of the 
pion photoproduction amplitude with the TriV phase-shifts 
for energies up to the first open channel. For energies 
above this opening, this theorem doesn't hold anymore. 
The energy region of interest here is in the second reso- 
nance energy region where the ttttN and rjN channels are 
opened Q- 

In order to understand the origin of this problem, we 
consider the pole approximation from Eq. ([39]) for ttN —^ 
ttN. In Fig.[7l the full solution of the Bethe-Salpeter equa- 
tion ([2]), using the simplified model from Ref. [5], is shown 
in red solid lines. The pole approximation from Eq. ([39]) 
is shown as the red dashed lines. Indeed, the resonance 
shape is well described by the residue term, while there 
are higher order terms ao, ai, ■ ■ ■ in the expansion of the 
full amplitude, which provide an almost energy indepen- 
dent background. 

The ttA'^ coupling is extracted at the pole position ac- 
cording to Eq. (|39p and given in the isospin basis (baryon 
first coupling convention) by g^NiSn) = + y/ 2 / 3\tt~^ n) + 
^/TJ3\tt^p) = 0.68 -f- 0.39 1. This means that the residue, 
expressed as a_i = g'^pj = jg^Arp exp(2 i</)7rAf) has an an- 
gle of 2 (pT^M = 60° which is significantly different from 
zero and causes the distorted shape of the A^*(1535) in 
Fig. [7] as compared to the classical resonance shape. 

In particular, the real part of the 6*11 amplitude (Re 
S'li) has a rise at energies around z = 1600 MeV, which 
is rather the shoulder of the A^*(1650) than part of the 
Af*(1535). In other words, the original fit from Ref. [S] tries 
to reproduce the narrow "valley" at Re z = 1550 MeV, 
while this structure is tied to resonance interference 1221 



Strictly speaking, the opening of the Ti~p channel which is 
just about 5 MeV above the tt'^p threshold already invalidates 
the applicability of the Watson's theorem. 
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Fig. 7. Simplified version without -kuN / form factors of the 
original model from Ref. [5] in -kN — ^ nN (red solid line) . The 
full model from Ref. [5] shows a slightly better agreement with 
the ttN PWA. The red dashed line shows the pole approxima- 
tion from Eq. (|39p . The data points are from the PWA of Ref. 

m. 



which is not included in the model of Ref. [5] . As a con- 
sequence, such a narrow structure requires a A^*(1535) 
pole close to the physical axis; indeed, a width oi F = 
— 2Inizo ~ 90 MeV has been found in Ref. [5] which is, 
although seen in some experiments [17| , rather at the lower 
limit of realistic values for the width [36| . 

Once wc consider £'□+ in pion photoproduction instead 
of partial waves in ttN scattering, the phases change, as we 
will see in the following. In photoproduction, the residue 
oLi is given by 



1 = 9-1 93^ 

6 



9i 



(40) 



with from Eq. The first line in Eq. ^ indicates 

that the residue can be decomposed in an effective photon 
coupling and the strong coupling gj from Eq. p9p for 
the final meson-baryon state in channel j. In the absence 
of genuine resonance states, the phototransition coupling 
g^ is given by the sum over all photon loops, including the 
strong transition strength g^ , as the second line in Eq. (|40p 
indicates. The sum is over the six channels in the particle 
basis as quoted in the first column of Table [T] The pole 
approximation in photoproduction is then given by the 
third hue in Eq. 

The sum of the photon loops is over the coupled chan- 
nels -kN , rjN, KA, and KS. This means the phototran- 
sition coupling is sensitive to all the strong couplings 
gi] they all appear in the transition, weighted by the re- 
spective photon loops. In particular, the residue a^i will 
obtain a new phase that is, in principle, very different from 
the phase of the residue a_i of TriV scattering. Thus, the 



phase of a resonance in £^0+ can reveal valuable informa- 
tion about the transition strengths into the different cou- 
pled channels and their relative phases and magnitudes as 
they appear in the sum in Eq. (|40p . 

As a consequence of its dynamical generation, the 
N*{1535) has large couplings to KA and KE [5], and it 
shows this feature in different theoretical approaches [l][5] . 
While the corresponding couplings cannot be accessed in 
direct experiments, they contribute to g-y, and thus, to the 
phase of i?o+. Studying Eo+, therefore, provides a valuable 
tool to check the magnitudes and phases of the strong 
couphngs to the dynamically generated TV* (1535) pole and 
thereby helps confirm or rule out models of dynamical 
generation. 

In contrast to i?o+, observables like cross sections do 
not provide any information of the phase, or only indi- 
rectly through interference with other partial waves. Study- 
ing Eo+ provides, thus, a much more sensitive test in pho- 
tonuclear reactions than in previous studies [23j . 

We can determine the phase of the residue al ^ similar 
to the case of ttA^ scattering. In the case of Eo+ , the phase 
is defined in the following way: we adopt a phase conven- 
tion in which a phase of zero degrees corresponds to a 
classical resonance shape in pi?o+ and „i?o+- This means 
a single maximum in Im£'o+ at the resonance position, 
while for KeEo+ a maximum below the resonance and a 
minimum above the resonance. Such a choice is fulfilled 
by the following definitions: 



pall = (-)■ 



(-)■ 



-iM 



N 



An z 



^aIi(n7r+) + iaIi(p^0) 



-iM 



N 



Att z 



V2 ^ 



^Jp.n) = arctan(lmp^„i. 



all iP'^ 



1 



(41) 



These definitions of pal^ and nO-li take into account the 
connection between Mf^ii and £'0+ as given by Eqs. (|36l38p . 
In Eq. (|41]) there is an additional minus sign that takes 
into account that a classical resonance shape is given by 
~l/{z — zo) and not l/{z — zo) [cf. Eq. ([40]) ]. In the calcula- 
tion of the phases as defined in Eq. (|4T|) an approximation 
is made by evaluating the d functions appearing in g-y from 
Eq. (|40|) at z = Rezo instead at the pole position itself, 

Z = Zq. 

The results for the phases of the A^*(1535) in photo- 
production, predicted by the model of Refs. [51123]. are 



i2 = 178°, ^7. = -9". 



(42) 



The phase for „i?o+ is —9"; indeed, this value close to 
zero reflects the classical resonance shape as observed for 
nEo+ (red solid lines) in Fig. [6l For pEo+ the phase is 
-t-lTS*^, and we observe an inverted resonance shape (red 
sohd lines) in Fig. [6l 

Thus, for the A^*(1535) with 7p or 771 initial state, 
there is a relative phase of nearly 180*^; this is in agree- 
ment with the findings of Ref. [23] where it was found 
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that the hclicity amplitudes A^^i^ and A^"/'^ have oppo- 
site sign. In Ref. [55] this has been interpreted as a suc- 
cess of the model of Ref. [5], because this opposite sign 
(or nearly 180*^ relative phase) is in agreement with the 
PDG values [35]. Yet, we have seen here that there is more 
than a relative sign - even with a 180° relative phase, the 
predicted individual multipoles pi?o+ and „i?o+ strongly 
deviate from the partial wave analyses as Fig. [S] shows. 

Summarizing the findings of this section, the predic- 
tion for i?o+ using the original model of Ref. [5] is seri- 
ously off the data. We have seen that there is already a 
potential problem in the description of the TriV ttN 
scattering data: as the interference of the iV*(1535) with 
the iV*(1650) is neglected, the model from Ref. [5] is forced 
to produce very small widths for the TV* (1535). Second, 
for i?o+, the effective photon coupling to the iV*(1535), 
which is a model prediction, produces a phase of the 
iV*(1535) in photoproduction which is in disagreement 
with the partial wave analyses of the Eo+ multipole. 

As discussed in SecjU we will reconsider the model for 
the A''* (1535) in Secs. l3.2[|3.3l bv allowing for the necessary 
degrees of freedom to solve the problems found in this 
section. This means the inclusion of the iV*(1650) as a 
genuine resonance, and a second genuine resonance that 
is allowed to replace the dynamically generated one if this 
is required by the fit. 

3.1.1 Factorization of the phototransition 

Discussed as case (IV) in Sec. 12.31 a further simplification 
of the phototransition amplitude can be carried out. We 
do not apply these further approximations in this study 
but discuss its consequences for the resonance phase. Case 
(IV) means an on-shell factorization of the photon loop. 
Then, in the absence of genuine resonances, the photo- 
transition amplitude, as tested in this subsection, can be 
written as 

AfflP, = (l-y^PG)-\% (43) 

where v!^ is given by the set of tree-level diagrams -fN — > 
MB with the meson-baryon pair MB on-shell. The pho- 
ton loop factorizes then into the form v!^ G with G being 
the factorized meson-baryon propagator. Eq. (|43p implies 
that the regularization of this factorized G is chosen to be 
the same as the G in the hadronic part, for each channel 
TriV, rjN, KA, KU. Such a choice, together with the factor- 
ization, has been realized e.g. in Ref. [5]. Additionally, an 
effective range expansion of the phototransition amplitude 
w/jj has been carried out in Ref. [2]. We test this simpli- 
fied amplitude within the model of Ref. [S] . Note that the 
hadronic interaction and the regularization (form factors) 
in the model of Ref. [5] are different. 

For the rjN, KA, and KS channels, the effective range 
expansion of vt!^ may be a sufficiently good approximation 
in the TV* (1535) region due to the proximity of the re- 
spective thresholds. However, the ttN threshold is farther 
away from the iV*(1535) region. While the effective range 
expansion is a good approximation up to 200 or 300 MeV 



above the ttN threshold, it is off by a factor of two in the 
A^*(1535) region. We have seen in the previous section [cf. 
Eq. (00])], that the resonance phase in photoproduction 
results from a subtle interference of the photon loops in 
the coupled channels; thus, this factor of two can cause 
large distortions in the iV*(1535) phase. 

Next, the full £^0+ amplitude is evaluated, using the 
effective range expansion for the on-shell factorized photo- 
transition, and the hadronic amplitude from Ref. [5j. With 
this treatment of the photon loops, the resonance angle is 
4>J = 33° which is more than 90° different from the value 
of Eq. dH]) of 178". In particular, the phase of 33° is not 
too different from the phase in ttN ttN. However, as 
has been shown here, this should not be considered as a 
solution to the phase problem (cf. Fig.lH]), but rather as a 
consequence of an oversimplified treatment of the photon 
loops. 

For the non-factorizcd photon loops used in this study, 
approximations have been made as well; however, as poin- 
ted out at the end of Sec. 12.4.11 these approximations 
only change the phase by a few degrees, at least within 
the present framework [cf. end of Sec. 12.4.1] . As a further 
test, in Ref. [55] the photon loop has also been calculated 
relativistically. In this case, the phase of the phototran- 
sition amplitude changes by 17° [13 which is still small, 
compared to the phase problem found in the previous sec- 
tion. 

3.2 Pion production at low energies (Fit 1) 

In the following sections, we present the fit results using 
the model of the present work introduced previously (see 
Sec. 12.21 for the hadronic part and Sec. 12.41 for the pho- 
totransition amplitude). Before considering the iV*(1535) 
and iV*(1650) region, it is instructive to study the low 
energy region. The low energy physics should not depend 
on resonances in the second resonance region, as required 
by chiral symmetry; this has been ensured by the use of 
derivative couplings for the bare ttNN* and 'jNN* ver- 
tices [cf. Eqs. ^ and ((57|) ]. Indeed, the corresponding 
parameters are very insensitive to the low energy region, 
as the fit shows. Thus, we have removed all contributions 
from the genuine resonances for the low energy fit. The 
only free parameters are then given by the four subtrac- 
tion constants (aTrjv, OjjJV, o-ka and axs) of the six coupled 
channels. Furthermore, those subtraction constants corre- 
sponding to the heavier channels, KA and KE, are very 
insensitive to the low energy region, as expected. We have 
fitted the S"— wave pion production, induced by photons 
and pions, up to 2; ~ 1.4 GeV in energy. The results for 
the subtraction constants are shown in Table [3] the re- 
sulting amplitudes are shown in Figs.lHldni andfTT] We 
refer to this low energy fit as "Fit 1" in the following. 

In Table [31 the parentheses indicate the weak sensi- 
tivity of the subtraction constants of the heavy channels. 
The j]N subtraction constant only has an impact at the 
higher border of the considered energy interval. The value 
of flTrJV = 2.65 is comparable to that from the original 
model of Ref. [5] of a„N = 2.0. 
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Table 3. Parameters of Fit 1 (low energy vrA^ nN and 
7A'^ — > ttN). The parentheses indicate less influential parame- 
ters. 



aKS (-3.80) aKA (3.80) a^N 2.65 a^jy 0.49 



PRL 87, 232501 (2001) 
PRC 55, 2016 (1997) [Im: gray band] 
Maid 2007 [EPIA 34,69] 
ChPT IZPhys C70, 483 (1996)] 
DR [PLB 399,13] .j^ 
present model 
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Fig. 8. (Fit 1, red lines); tt" photoproduction close to thresh- 
old. The n'^n threshold is indicated with an arrow. Experimen- 
tal analyses: black (gray) data points from Refs. [JS] (|50)). 
Gray band from Ref. [50]. Theory: Blue double dashed dotted 
lines: MAID 2007 [30], brown dashed lines: Ref. [51]. Dashed 
dotted lines: ChPT calculation from Ref. [52] (see also the more 
recent work of Ref. [53] and the analysis of Ref. [56] ) . 



The 7r°p photoproduction close to the ttN threshold of- 
fers a very sensitive test of the phototransition amplitude. 
Precise data exist for i?o+(7r°p), that include the cusp ef- 
fect from different pion and nucleon masses. As the present 
model is formulated in the particle basis rather than the 
isospin basis, these effects can be taken into account in the 
present work. The result of Fit 1 for £'o+ (tt'^p) is shown in 
Fig. [5] with the red solid lines. The experimental analyses 
are from Ref. [50] . Note that for Im Eq+ , only the experi- 
mental band from Ref. [SD] is shown. In the other analysis 
from Ref. [32], a comparable value for Im Eq+ is obtained. 
In Fig. [8l also results from MAID are shown |30n5l] and 
from a ChPT calculation [5^ . 

The first thing to note is that i?o+(7r*'p) is almost in- 
dependent of the fit parameters. This is because the tree- 
level contribution and the one-loop contribution (i.e., the 
photon loop without rescattering) do not depend at all 
on the subtraction constants; the two-loop amplitude is 
orders of magnitudes smaller so close to threshold. 



Then, Im i?o+(7r*'p) is entirely given by the ir^n one- 
loop amplitude. The strongly energy dependent cusp struc- 
ture in Re Eq^ (t^^p) is entirely given by the dispersive part 
of that loop. The tree-level photoproduction diagrams as 
well as the loop contribute to an almost energy indepen- 
dent background in Re Eq+Itt^p). 

As Fig. [5] shows, the photon loop evaluated in Sec. 
12.4.11 indeed predicts the correct energy dependence for 
both real and imaginary parts of Eq^ (t^'^p) ■ This is a good 
test that our phototransition amplitude provides a realis- 
tic picture close to threshold. 

We had to shift the theoretical result for Re Eq+ (tt^p) 
by —0.59 X 10~^m~+, as indicated in the figure with an ar- 
row. To judge the size of this amount, we consider the dif- 
ferent contributions to Re Eo+{tt'^p). From the tree-level 
diagrams of Fig. [3] only the direct and crossed nucleon ex- 
change (b) and (d) contribute. For ^-wave photoproduc- 
tion, these diagrams are typically one order of magnitude 
smaller than the KroU-Ruderman term (c), which vanishes 
identically for neutral pion photoproduction. This already 
shows the sensitivity of £^0+ to higher order correc- 
tions. 

Still, these subleading contributions from Fig. [3] (b), 
(d) are five times larger (—2.47 x 10~^m~+ ) than the ex- 
perimental value at the 7r+n threshold of ~— 0.5x 10^'^to~+ . 
The photon loop contributes with 4-2.5 x 10~^m~|; this al- 
most fully cancels the tree-level contribution, so that the 
final theoretical result is around 0.59 x 10~'^m~+ larger 
than the experimental value. 

Thus, even the subleading tree-level diagrams (b), (d) 
vanish through cancellation by the photon loop. A part 
of the remaining discrepancy of 0.59 x 10~^m~l comes 
from the slightly inconsistent treatment of the tree level 
diagrams and the loop contributions, since the tree level 
is treated fully relativistically, while for the loop a non- 
relativistic framework is used. 

Neutral pion production close to threshold has been 
calculated in the framework of chiral perturbation the- 
ory [51[53[53Mi|. At next-to-leading (NLO) order, both 
the anomalous magnetic moment coming from the tree 
level diagrams and the triangle diagram contribute. The 
nonanlytical piece from the triangle diagram at NLO is 
given by 



AEo+{7t"p) 



128 7rVr 



(44) 



In Fig.[9l this pion mass dependence (dashed line) is com- 
pared to the contribution from the triangle diagram in the 
present formulation (solid line), given by Fig. Hj^a). For 
this comparison, we have switched off the final state in- 
teraction MB MB, i.e., the shaded circle in Fig.[3I^a) is 
given by the on-shell factorized Weinberg- Tomozawa term. 
Second, we consider only the ttA^ loop, and the tiny con- 
tribution from the photon coupling to intermediate KY 
states is switched off. 

The present result has the same /i (= m^^ /Mj^) de- 
pendence as the result from Eq. (j^ (it does, though, 
include also higher order pieces) and both results seem to 
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Fig. 9. Pion mass dependence of the triangle diagram. Solid 
line: Contribution to Eo+{-K^p) from the triangle diagram of 
Fig. |4ja) [meson pole term]. Dashed line: Nonanalytic piece of 
the NLO contribution given in Eq. (|44|) . 



be in reasonable agreement as Fig. [5] shows. Yet, compar- 
ing the result at the physical pion mass (/i = 0.15), what 
seems to be a small difference between both results, of 
0.3 X 10~^m~+. is of similar size as the previously stated 



discrepancy of 0.59 x 10 



-3, 



shown in Fig. [5] We can, 



thus, conclude that although the present calculation can- 
not precisely describe the data of Fig. [U those discrep- 
ancies are much smaller than the NLO contribution in 
the chiral expansion 3 x 10~^m~+); neutral pion pro- 
duction close to threshold is difficult to describe as it is 
sensitive to higher-order corrections. 

As final remark on i?o+(7r*'p) near threshold, we noted 
that in the present work the real part of the photon loop 
is finite and fixed by gauge invariance, and there is no 
freedom from a regulator to adjust the model to the data. 
In principle, the discrepancy in i?o+ discussed above may 
be eliminated by a proper adjustment to the data of the 
transverse contact current T'' in Eq. ([8|), which has been 
set to zero in the present work. However, we have cho- 
sen not to do so in order to avoid introducing more free 
parameters in the current model. 

The outcome of Fit 1 for -kN ttN and jN — > ttN, 
fitted both in and 6*31, is shown in Figs. [TOl and [TT] 
The results show a good global agreement with the "data" 
from the partial wave analyses, given that uka and axs 
are almost insensitive to this energy region; the fit thus 
describes, essentially with the two free parameters a^^r 
and QriNi ten different data sets. Still, there are remaining 
discrepancies in Re Su for ttN — > ttN close to the wN 
threshold, and for Im £^0+ on the proton. It is clear that 
the present model cannot deliver higher accuracy results 
at this point, because the interaction has been limited to 
the lowest order chiral interaction. It is known (see e.g. 
Ref. [59]) that close to the ttN threshold, higher order 
terms in the chiral interaction are needed to deliver a con- 
sistent threshold behavior within the present scheme of 
hadronic interaction. As this problem is well-known and 
our interest is different in the present study, we do not try 
to solve these discrepancies and refrain from calculating 
the scattering lengths. 
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Fig. 10. (Fit 1) Low energy pion scattering. Analysis without 
resonances. Data from partial wave analyses as in Fig. [T] 
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Fig. 11. (Fit 1) Low energy pion photoproduction. Analysis 
without resonances. Data from partial wave analyses as in Fig. 
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In the next section, we will include higher energies in 
the fit. It is then not possible to fully maintain the qual- 
ity of Fit 1 for the low energy region. We have already 
seen that one reason is the importance of higher orders in 
the chiral interaction. Second, the on-shell factorization 
scheme [cf. Sec. 12. 1| leads to the appearance of subtrac- 
tion constants, which are independent in energy. This en- 
ergy independence is, of course, only an approximation, 
because higher order interactions induce more divergent 
loops, which require a multiple subtraction, i.e. a polyno- 
mial in energy instead of a constant. We thus have to rec- 
ognize that the present scheme can only fit the data within 
an energy window, in which the subtraction polynomials 
are locally given by a constant value to a good approxi- 
mation. Yet, while the model is limited at this point, one 
can stiU fit the entire A^*(1535) and iV*(1650) region to 
a sufficient precision while maintaining at the same time 
the main features of the low energy behavior of Fit 1. 



3.3 Pion production in the Second Resonance Region 
(Fit 2) 

In this section, the main results of the present work are 
presented. After exploring the low energy region in the 
previous section, we now turn to the second resonance re- 
gion employing the full model including the genuine res- 
onances. The fitted data are the 5*11 amplitude in ttN —s- 
ttN plus the S31 amplitude below the region of the Z\(1620) 
(this resonance is not included in the present model). Si- 
multaneously, the 5*1 1 partial wave analysis data for _Bo+ 
on the proton and on the neutron plus the 6*31 Eo+ mul- 
tipole below the zi(1620) region is included. For all data 
on Eo+ , we have only included the imaginary part of the 
amplitude; Re £'0+ is in all cases calculated but not in- 
cluded in the fit. This is because in the previous section we 
have seen that there is some theoretical uncertainty from 
higher order tree-level diagrams, which all contribute to 
the real part (up to unitarity corrections). In the fit, we 
have given more weight to the second resonance region, 
but also some weight to the low energy region to maintain 
the main features of the low energy Fit 1 discussed in the 
previous section. 

The parameters of the solution are shown in Table ID 
There are four subtraction constants, and for each of the 
two genuine resonances, we have four bare hadronic cou- 
plings, two bare electromagnetic couplings, and one bare 
mass. We refer to this fit as "Fit 2" in the following. 

The resulting amplitudes are shown Figs. [T2l and [T3l 
At low energies, the results are slightly worse than for 
Fit 1 (compare to Figs. [TOl and [TT|) . while the threshold 
region for i?o+(7r''p), shown in Fig. [8] for Fit 1, is nearly 
unchanged. This is because the influence of resonances 
vanishes at the t:N threshold, as required by chiral sym- 
metry. Second, the subtraction constants appear first at 
the two-loop level in photoproduction, which is very small 
close to threshold. 

As already discussed at the end of the previous section, 
one cannot maintain the quality of Fit 1 , once the second 
resonance region is included in the fit. The reasons have 
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Fig. 12. Data points: SES irN ttN partial wave analysis 
from Ref. [45] [FA07]. Red line: Joint analysis of ttN nN 
and jN -> tvN (Fit 2). 




-40 i-L 



1200 1400 1600 
z [MeV] 



1200 1400 1600 
z [MeV] 



Fig. 13. Red lines: Joint analysis of ttN nN and 7A'' 
(Fit 2). Data points and other lines: Labeling as in Fig. 
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Table 4. Parameters of Fit 2 (TriV TriV and 77V -> TriV, 
entire energy region). Subtraction constants, bare strong and 
electromagnetic couplings, and bare masses [cf. Eqs. (|3l I27[l]. 
The parentheses indicate less influential parameters. 
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been pointed out: higher orders in the chiral expansion 
are not included in the present scheme; closely connected 
to this, the subtraction polynomial is a constant in en- 
ergy in the present model which restricts the freedom of 
the fit; third, the ttttA^ channel is not included. However, 
while higher order corrections will necessarily deliver a 
more precise fit, the qualitative results from the current 
model should not change. In this study, we are interested 
in the interplay of genuine and dynamically generated res- 
onances, and the precision of the current model is sufficient 
for this discussion. 

In any case, in the second resonance region. Fit 2 de- 
livers a fair data description while maintaining the main 
features of Fit 1 at low energies. The results shown in 
Figs . fT2l and fT3l are much better than those from the orig- 
inal model of Ref. [5]; in particular, the phase problem, 
pointed out in Sec. l3.1l has been solved. Note also that the 
results are in fair agreement with Re Eo+, although only 
Im £^0+ has been included in the fit. The present results 
are in better agreement with data in the -/V*(1535) and 
A^*(1650) region than in the previous work [5] within the 
framework of UxPT. In the S31 partial-wave state, there 
are, of course, deviations in the region of the Z\(1620) as 
this resonance has not been included. 

Note the appearance of a cusp in Re 5*11 in ttN irN, 
while the imaginary part docs not show such a pronounced 
cusp. In photoproduction, this is different: there are strong 
cusps both in the imaginary and real parts of i?o_|_ in 
Sii{p) and S'ii(n) as Fig. [T51 shows. This different func- 
tional behavior is allowed as Watson's theorem does no 
longer hold above the ttttA^ threshold. Fit 2 reproduces 
these different functional forms. In fact, the previous SES 
analysis of i?o+ from Ref. [29] (gray data in Fig. [T3|) did 
not clearly show the sharp cusp in Im Eq^ Sii{p), which 
first led to major concerns, because in the present model 
of dynamical generation of the 7V*(1535), that cusp is un- 
avoidable and it is a very stable feature. However, the 
most recent analysis from Ref. [45] (black data in Fig. [T3)) 
clearly reveals, that this cusp is indeed present, and very 
pronounced. This is shown in greater detail in Fig. 1141 

In Ref. [60] the role of the rjN cusp in TriV — > TriV 
has been discussed. The iV*(1535) could be dynamically 
generated but at the cost of a worse data description 
and the appearance of a sharp pronounced cusp. It has 
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Fig. 14. Detail of Fig. [13] in the iV*(1535) region. Note the 
form of the cusp at the T^iV threshold is quite different from 
the TriV TvN reaction. 

then been argued that a genuine iV*(1535) is preferred 
due to the better description of the 5*11 partial wave. 
While in the present study, sharp cusps indeed appear, 
the data description is satisfactory; the dynamically gen- 
erated iV*(1535) is characterized by a true pole as will 
be seen in Sec. H] and the cusp will be shown to be an 
interplay of physical and hidden poles, which necessarily 
appear. Note also, that in the latest version of the Jiilich 
model |20] a sharp 77iV cusp appears despite of the fact that 
the iV*(1535) is introduced as a genuine resonance state 
in that model. The worse data description in Ref. [60] . 
for the case of the dynamical generation of the iV*(1535), 
may be tied to the missing KA and KS channels in that 
model; including these channels in the Jiilich model may 
help settle the issue of the nature of the iV*(1535). 

As argued in Sec. 12.4.11 wc have used only the trans- 
verse part of the meson pole term for the phototransi- 
tion loop amplitude. In contrast, the contribution from the 
baryon pole term [cf. Fig. [J has been neglected. Indeed, 
we have checked that the contributions from this term are 
small. They arc smaller than 10 % for the photoproduc- 
tion amplitude on the proton, and slightly larger than 10% 
for the photoproduction amplitude on the neutron. In any 
case, one can safely neglect these contributions at the level 
of precision we arc working in this study. 

While both the iV*(1535) and the iV*(1650) are well 
reproduced by Fit 2, both in photon- and pion- induced 
reactions, we cannot conclude from the amplitude on the 
real z axis, which role the genuine poles play. The val- 
ues from Table m already indicate that the second genuine 
resonance has probably moved far into the complex plane 
(large bare couplings and mass). After all, only a study of 
the complex z plane on unphysical sheets can resolve this 
issue. This will be carried out in the next section. 

Here, we anticipate some results of the next section. 
Indeed, in our final results (Fit 2), the iV*(1535) is dy- 
namically generated, while one of the genuine resonances 
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is responsible for the A^*(1650) and the other one provides 
a background that varies very slowly with energy. This 
"background pole" is far in the complex plane. As a test, 
we have refitted the data without the genuine resonance 
leading to this pole, but the becomes worse, indicating 
that the remaining free parameters of the model cannot 
compensate for the absence of the background pole. In 
order to see in which reaction this background is more 
important, we have performed further tests without sec- 
ond genuine resonance, fitting to (a) only ttN ttN and 
(b) only jN ttN. In case (a), we can obtain a good 
fit, with a dynamically generated iV*(1535) and a genuine 
TV* (1650), without the need for an additional background. 
In case (b), however, we cannot obtain a good fit; in the 
best fit, the only remaining genuine state even becomes a 
background pole (i.e. moves far into the complex plane), 
instead of describing the N*{1650). 

Thus, we can conclude that, particularly in photopro- 
duction, an additional background is required on top of 
the N*{1535) and TV* (1650) and the background gener- 
ated by the rescattering provided by Eq. ([2]) from V^^ . 
Such an additional background in photoproduction could 
be provided by diagrams which are not explicitly included 
in the present study, like vector meson ^-channel exchange 
with anomalous photon couplings. 

Finally, we have also studied the importance of the KA 
and KE channels for the formation of the iV*(1535). It 
is known that these channels are responsible for the nec- 
essary attraction to generate the TV* (1535) [UOE]. In a 
simple model without genuine resonances, we have fitted 
the ttN —f ttN reaction, only in Sn and only around the 
resonance position. Additionally, we have increased fx in 
Eq. ([T|) from its original value of /k ~ 1-22 /^r = 113.5 
MeV (a more precise value of fx = 1-193 has been 
reported recently in Ref. [35]); this models a weaker cou- 
pling to the strangeness channels KA and KE. Within a 
reasonable range for the subtraction constants (maximum 
3 to 4), we could get a resonant shape of the iV*(1535) 
for values of Jk up to 150 or 160 MeV. For larger Jk, 
the resonance fades away. A maximum value of fx = 150 
MeV has been regarded in Ref. [40] as a reasonable limit. 

Thus, we can conclude that the dynamical generation 
of the TV* (1535) indeed requires a sufficiently strong cou- 
pling to the KA and KE channels. Since in the present 
approach, the coupling strengths to these channels are 
given by the SU(3) Lagrangian which are strong, the dy- 
namically generated A^*(1535) appears as quite a stable 
structure in the various refits discussed in this work. 

More data can be included in the fit that allow to fur- 
ther test the SU(3) structure, such as those from the pion- 
and photon-induced rjN, KA, and KE production, as well 
as from the corresponding electroproduction processes. As 
we have 18 free parameters altogether, an inclusion of 
those data is appropriate to impose further constraints 
on the model. The results for these observables will be 
presented in Ref. [HI] ; the current Fit 2 already delivers a 
good qualitative agreement for these observables. 



4 Discussion 

4.1 Pole positions and residues 

Poles and zeros of the amplitude in the complex plane 
of the scattering energy z = s^^^ determine the global 
appearance of the amplitude on the physical axis. The 
hadronic amplitude can be analytically continued to the 
complex 2;-plane. There are two different Riemann sheets 
for each channel ttN, r]N, KA, and KE [22] . 

The first sheet of the propagator G is defined by evalu- 
ating G for complex z, while the second sheet is obtained 
by adding twice the discontinuity of G along the right- 
hand cut, 



GW(2) = G(z), 

G(^)(z) = G(.) + 2^ 
47r z 



where M is the baryon mass of a given channel and 



> 

on 



Qon 




if Im Qo 
else 



< 



(45) 



(46) 



is the on-shell relative momentum. Eq. pS]) ensures that 

has the cut along the right-hand side. With this pre- 
scription, both G^^-* and G'-^-' have the cut along the posi- 
tive physical axis and are analytically connected with each 
other along these cuts. 

The various sheets of the scattering amplitude T are 
induced by the replacement of G from Eq. ([2]) with G*-^-* or 
G^^^ . This amounts, for the four meson-baryon channels, 
to 16 sheets of the scattering amplitude (in this section, 
we work in the isospin limit, so that the 6 channels in 
the particle basis reduce to the four channels ttN, i]N, 
KA, and KE). Out of these 16 sheets, only a few are di- 
rectly connected to the physical axis. Directly connected in 
this sense means connected without having to turn around 
branch points to reach the physical axis. This issue is dis- 
cussed in detail in Ref. [22j . 

In the following, the physical sheet of T induced by 
G^^^ is labeled 1, the second, or unphysical sheet, induced 
by G(2), is labeled 2, for a given channel. For example, 
in the channel ordering ttN, rjN, KA, KE, sheet 1111 
is the physical sheet with respect to all channels. It is 
free of poles. Sheet 2211 is the sheet that is given by the 
unphysical sheet of wN, rjN, and the first sheet of KA and 
KE. In the following, we concentrate on the lower z half 
plane. The properties of the amplitude in the upper z half 
plane are analogous to those of the lower half plane and 
T*{z) ~ T{z*) (Schwartz's refiection principle). 

As discussed in Ref. [22], for a given channel, the first 
sheet in the lower z half plane is directly connected to the 
physical axis below the threshold z < mi + Mi. In con- 
trast, the physical axis above threshold z > rrii + Mi is 
directly connected to the second sheet in the lower z half 
plane. For example, the sheet 2111 is directly connected 
to the physical axis for -I- < z < + M^. Thus, 
there are four combinations of sheets that are directly con- 
nected to the physical axis, 2111, 2211, 2221, 2222. The 
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Table 5. Positions and couplings g^-p [cf. Eq. (|39p ] of poles 
on different sheets. VS means the virtual state below the irN 
threshold, BG means the background pole. The physical poles 
are underlined [cf. discussion in Sec. 14.2) . The virtual state 
and the A*'* (1535) appear dynamically generated, while the 
A'^*(1650) and the background pole result from genuine reso- 
nance states. 





Position [MeV] 


n 

bfTT p 


Sheet 2111 








VS 


1031 - 


203 i 


-0.51 + 1.58 i 


Af*(1535) 


1647- 


103 i 


-1.55 + 1.40i 


AT* (1650) 


1872 - 


57 i 


0.91 + 2.64i 


Sheet 2211 








iV*(1535) 


1608 - 


175 i 


3.35 + 1.82 i 


AT* (1650) 


1645 - 


105 i 


-1.83 + 1.88i 


BG 


1545 - 


545 i 


-0.78 + 3.52 i 


Sheet 2221 








A^*(1535) 


1538 - 


139 i 


1.42 + 0.46 i 


Af*(1650) 


1655 - 


59 i 


-0.89 + 0.48 i 


BG 


1837- 


800 i 


0.31 + 2.39 i 


Sheet 2222 








A^*(1535) 


no pole 




iV*(1650) 


1662 - 


72 i 


-1.03 + 0.12i 


BG 


2129 - 


1289 i 


0.33 + 2.26 i 
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corresponding pieces of the physical axis for these sheets 
are indicated in Fig. [15] with the bold red lines. 

Usually, poles are only searched on one combination of 
the sheets [5111 3| using the prescription for the propagator 
G, 



G,{z) 



G?\z) 



if Re z < rrii 
if Re z > nii 



M,. 



(47) 1 



for the channels i with meson and baryon mass nii and 
Mi, respectively. The prescription given by Eq. (|47)) _is a 
sensible one, because the Riemann sheet induced by G{z) 
is the closest one to the physical axis. Poles on other sheets 
are expected to have much less impact on the amplitude 
on the physical axis. 

However, studying all four combinations 2111 to 2222 
instead of using G from Eq. (|47p reveals also virtual states 
and hidden poles, that are not found using the prescription 
from Eq. (|47p . Those sub- and above-threshold resonances 
can have a large influence on the physical amplitude as wc 
will see in the following. 

In Fig. I15[ we show contour plots of the amplitude for 
the contours Re T{z) = (solid lines) and Im T{z) = 
(dashed lines). It is easy to see that those lines intersect 
at poles (red circles) and zeros (blue crosses) of the am- 
plitude. In the following, we refer to this kind of represen- 
tation of a complex function as "GauB plot" . 

The pole positions and the couplings to the Tr~p chan- 
nel, defined in Eq. ([39]) . are listed in Table[5l The coupling 
strengths to isospin 7 = 1/2 can be obtained through 
91=1/2 = \f^T^9-K-p up to small isospin breaking from 
different masses which is only important for the cusp ef- 
fect in _Bo+ {''^^V) 7 shown in Fig. [51 The couplings to the 
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Fig. 15. "Gaul3 plot" of the Riemann sheets 2111 to 2222. The 
contours Re T = (solid lines) and Im T = (dashed lines) 
intersect at poles (red circles) and zeros (blue crosses) of the 
amplitude. The part of the physical axis directly connected to 
the respective sheet is indicated in bold red. 
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other channels gr/N^ gKA, and gKS are not listed in Table 
m although they have been calculated. They will be fully 
quoted in Ref. [61], once the observables in the ijN, KA, 
and KE channels are included in the fit. 

There are two genuine poles in the model. One has 
become the TV* (1650) in the solution, the other one has 
become a background pole that lies far in the complex 
plane and provides an almost energy independent back- 
ground. The virtual state VS (discussed below) and the 
A^*(1535) appear as dynamically generated. 

The couphngs of the 7V*(1535) on sheet 2211 to KA 
and KS are l^^-o^j = 4.3 and \gK+s- \ = 2.3. The large 
coupling to KS indicates that the A^*(1535) appears as 
a quasibound KS state, in qualitative agreement with 
the models of dynamical generation from Refs. [Tl[2] and 
[5]. Yet, the situation is more complicated here, because 
the interference with the A^*(1650) is nontrivial, as will 
be discussed in Sec. 14.21 Second, the coupling strengths g 
of the A^*(1535) on sheet 2211 -which is the pole to be 
compared with the pole of Ref. [5]- are larger than those 
from Ref. [S], for all channels. However, this is a simple 
consequence of the fact that the A^*(1535) is much wider 
here than in Ref. [5]. In order to achieve a comparable 
resonance shape on the physical axis at Im z = 0, a pole 
located farther in the complex plane needs to have a large 
residue according to Eq. (|39p . 



4.2 Discussion of four sheets 

In Fig.[T5l consider first sheet 2211, which is connected to 
the physical axis in the range + Mn < z < mx + Ma- 
This sheet shows the global behavior of the solution: the 
iV*(1535) is dynamically generated with a large width, the 
first genuine pole is identified with the A^*(1650), and the 
second genuine pole has moved far into the complex plane, 
providing a background that varies very slowly with en- 
ergy on the physical axis ("background pole"). The latter 
pole models additional background processes that are not 
explicitly included in the present model. Thus, there are 
several poles and zeros on sheet 2211. One of the zeros 
is situated in between the Af*(1535) and Af*(1650). This 
zero is also found in Refs. [35] and [52]. The other zero is 
near z = 2000 - i MeV. As it lies above the KA and KS 
thresholds, it is physically not observable; a zero on that 
position could only be observable on sheet 2222. See also 
a discussion in Ref. [20j on a zero on the physical axis in 
the ttN — !■ TjN transition in the Sn state. 

What are the physical implications of the poles on 
sheet 2211? Sheet 2211 is connected to the physical axis 
in the range m,, -I- Mm < z < tuk + Ma- The pole of the 
A^*(1535) is within this energy window. We can identify 
this pole with the physically observable N*{1535)- The 
iV*(1535) pole lies relatively far in the complex plane at 
Im z = —175 MeV; in the recent work of Ref. [SI], a 
similarly large imaginary part for the iV*(1535) pole has 
been found (Im zq = -191 MeV) while the iV*(1535) pole 
found in Ref. [22] is much closer to the physical axis (Im 
zo = -64.5 MeV). 



The real part of the pole position of the iV*(1650) on 
sheet 2211 is above rriK + Ma- Thus, the TV* (1650) on 
sheet 2211 is rather an above-threshold resonance. Only its 
low-energy tail is visible on the physical axis. The physical 
A^*(1650) has to be searched for on sheet 2221 instead (see 
below) . 

Sheet 2111 in Fig.[Tn]is connected to the physical axis 
in + Mn < z < -|- M^- Below the ttN threshold 
and far from the real axis into the complex plane, there is 
a pole on sheet 2111 as Fig. [T51 shows. We have found this 
virtual state in the Sn partial wave also in other models 
[22] : as discussed below, it is connected to the sharp rise 
of Re Sii at the ttN threshold and seems to be required 
by the partial wave. However, it is not clear if this state 
is genuine or a "forced" pole that mocks up the u- and t- 
channel subthreshold cuts that are not explicitly included 
in the present model |63| . 

The sheets 2221 and 2222 are also shown in Fig. [15l 
They are connected to the physical axis within the ranges 
rriK + Ma < z < rriK + Me and rriK + Ms < z < oo, 
respectively. The structure of sheet 2221 is similar to that 
of 2211. The pole of the 7V*(1650) on sheet 2221 is directly 
connected to the physical axis, for which we can identify 
it with the physical TV* (1650). The physical iV*(1535) and 
iV*(1650) on their respective sheets are highlighted in Ta- 
ble [5] As discussed before, the other TV* poles on other 
sheets always appear either as a sub-threshold or above- 
threshold resonance. Yet, also these secondary poles are 
important, because their tails can be visible on the phys- 
ical axis, and their interference with the physical poles is 
important. 

The structure of sheet 2222 is quite different from those 
of 2211 and 2221. Sheet 2222 is the sheet connected to the 
physical axis above the KS threshold. The A^*(1535) has 
disappeared on that sheet and the background pole has 
moved even farther into the complex plane (out of the 
plotted range, see Table [5]). 

The disappearance of the N*{1535) pole on sheet 2222 
is an interesting fact; the model of dynamical generation 
of the N*{1535) from Ref. [5] shows the same behavior. 
The absence of the N*{1535) at high energies has implica- 
tions for the concept of sub-threshold resonances that are 
sometimes used to fit reaction data: for example, in Ref. 
[64] . the role of the TV* (1535) in the reaction n~p (j)N 
was discussed. There, the A^*(1535) was used as a sub- 
threshold resonance to explain the (j>N cross section, i.e., 
the resonance was extrapolated several hundreds of MeV 
above its position. It was already argued in Ref. |65| , that 
one should rather use the full energy dependent meson- 
baryon amplitude instead of extrapolating the resonance; 
with the detailed study of the Riemann sheets done here, 
we can further sharpen this statement: the influence of the 
N*{1535) has just disappeared completely above the KS 
threshold. 

The disappearance of the A^*(1535) above the KS 
threshold is, of course, tied to the present model. Yet, 
even in cases the TV* (1535) does not disappear from the 
2222 sheet, due to its strong coupling to the strangeness 
channels (implicitly assumed when used to fit the pro- 
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ttN and -yN ttN 
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Fig. 16. Pole approximations Tp^ (piecewise defined lines) of 
the full solution for Re Sii. The pieces of the PA are limited 
by the thresholds of the channels nN, -qN, KA and KE. 

duction), its position will unavoidably change drastically 
on the 2222 sheet compared to its position on the 2211 
sheet, where it is observed. This is a model-independent 
behavior, and makes any kind of models using phenomcno- 
logical sub-threshold resonances questionable. 




1400 

Fig. 17. Cusp of the physical amplitude (schematically). 
Above the 77A'' threshold, the amplitude is dominated by the 
physical A''* (1535) pole on sheet 2211. Below the rjN thresh- 
old, the amplitude is dominated by the hidden A* (1535) pole 
on sheet 2111. 



4.3 The impact of poles on the physical axis 

The various poles on different sheets have different im- 
pact at the physical axis. Consider the pole approxima- 
tion from Eq. p9|) . that provides the leading term in the 
Laurent expansion around the pole position. In Fig. [TBI we 
shows the real part of the Sn amplitude together with the 
expansions from Eq. (|39p . 

Below the rjN threshold {z < 1487 MeV), two curves 
are shown in Fig. 1161 The dotted curve originates from the 
virtual state below the ttN threshold on sheet 2111 shown 
in Fig. [15] [cf. Table [S] . It almost saturates the amplitude 
close to the ttN threshold, which indeed shows that the 
low energy region is dominated by this virtual state. 

At the rjN threshold, the Re Su amplitude shows a 
characteristic cusp. The shape above the 7]N threshold is 
well approximated by the poles on the sheet 2211, while 
below the rjN threshold, it is well described by the "hid- 
den" A^*(1535) pole on the 2111 sheet, as indicated in Fig. 
[injwith the dashed line. Before studying the other pole ap- 
proximations above the rjN threshold, plotted in Fig. [TBI 
we further discuss the cusp. 

In Fig. [TT] the cusp and its origin are schematically dis- 
played. The physical amplitude at Im z = 0, around the 
cusp, is connected to two different Riemann sheets, 2111 
and 2211. The 7V*(1535) poles are at different positions 
on these sheets, and their influence on the amplitude, in- 
dicated with thick arrows, results in a cusp. This explains 
naturally the appearance of the cusp. 

The physical axis above the rjN threshold is divided 
into three pieces, separated by the KA and KU thresh- 
olds. The pieces are connected to sheets 2211, 2221, and 
2222, respectively. We have summed the contributions from 
the poles from Fig. [15] for each of these sheets according 



to 

J^S = E^ (48) 

where (?) indicates the amplitude on sheet i =2211, 2221, 
2222 and the sum is over the poles on a given sheet. The 
different ttN — > ttN residues and pole positions on a given 
sheet (z) are indicated as a-'_^ and Zq. It is important to 
sum over the different poles. In particular, the combina- 
tion of the Af*(1535) and A^*(1650) poles is important. 
This is because the A'^*(1535) provides a strongly energy- 
dependent background for the iV*(1650) and vice versa. 
In other words, the interference between these two res- 
onances is responsible for the shape of the partial wave 
amplitude; they cannot be treated separately. This issue 
has been extensively discussed in Ref. [22]. 

(i) 

The resulting pole approximations Tp^ are shown in 
Fig. [M] with the solid lines. They arc different for each re- 
gion of the physical axis between the channel thresholds. 
The results provide a good reproduction of the energy de- 
pendence of the amplitude, in particular the resonance 
shapes of the TV*(1535) and iV*(1650). Yet, for each piece 
of the physical axis above the r/N threshold, the pole ap- 
proximations are still off the full solution by backgrounds 
that slowly vary with energy. These almost constant off- 
sets, that are different for each piece, come from higher 
order terms in the Laurent expansion around the pole po- 
sitions, flo, fli according to Eq. (j39p . 

Finally, let us mention our results of the pole search 
in the model of Ref. [5]. In that reference, only the dy- 
namically generated pole of the iV*(1535) on sheet 2211 
has been found. However, we can perform the same de- 
tailed pole search as carried out for the present model. 
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Then, one indeed finds more poles in the model of Rcf. 
[5]. As in the present study, there is another pole of the 
A^*(1535) on sheet 2111. Similarly to the present model, 
the cusp structure at the t]N threshold can be well de- 
scribed by these two poles. Also, the virtual state below 
the ttN threshold is present in the model of Ref. [5]. In 
the same way as found here, that state is responsible for 
the sharp rise of the real part of the amplitude close 
to threshold. Furthermore, the model of [5] has another 
dynamically generated resonance on the sheet 2111 far in 
the complex plane at zq = 1657 — 267 i MeV that is re- 
sponsible for some of the structure along the physical axis 
between the irN and rjN thresholds. 

4.4 The phase problem revisited 

In light of the analysis of Fit 2, the phase problem in the 
model of Ref. [5] , found in the present study and discussed 
in Sec. 13.11 can be revisited. As we have seen in Figs. [121 
[T51 and [m the amplitude of Fit 2 shows no phase prob- 
lem any more on the physical axis. In the previous section, 
we have seen that the physical axis is dominated by the 
different iV*(1535) and iV*(1650) poles on the different 
sheets. In particular, one has to consider the interference 
of resonances, and individual contributions make no phys- 
ical sense. In view of this, one should be cautious to quote 
phase angles for the A^*(1535) alone. 

Yet, the interfering poles dominate the energy depen- 
dence in — > ttN and ttN ttN, as has been discussed. 
In particular, the phases of the resonances, together with 
the residue strengths and pole positions, lead to very dif- 
ferent resonance shapes as a comparison of Figs. [T2l and [T3l 
shows. Most noticeably, the different strengths of the rjN 
cusps in the real and imaginary parts of the 7 ttN and 
ttN — > ttN amplitudes can be explained by Fit 2: the in- 
terfering resonances with different phases modify the am- 
plitudes in such a way, that the different functional forms 
turn out naturally. 

In particular, the A* (1535) shape appears narrower in 
Eq^ than in nN — > ttN (see e.g. the narrow width found in 
Ref. [30]). However, in the present framework, this could 
be explained naturally by the photon coupling to the dy- 
namically generated A^*(1535) which induces a different 
phase [cf. Eq. (jlQ]) ] on the couphng constant. With such 
a new phase in photoproduction, together with the res- 
onance interference, the A^*(1535) can naturally appear 
narrower on the physical axis, while its pole position is, of 
course, still the same to that in ttN —* ttN. 

5 Conclusions 

The A^*(1535) has been previously described as a purely 
dynamically generated resonance from the unitarized low- 
est order chiral interaction in SU{3) coupled channel dy- 
namics. This concept has been tested in a variety of pion- 
and photon-induced reactions. In this study, we carry out 
a further test, which is more sensitive because it is directly 
tied to the amplitudes instead to cross sections. This is a 



very sensitive test to get further insight into the nature of 
resonances. 

The simultaneous study of the reactions ttN ttN 
and 7A —t nN for the 6*11 and S31 partial waves reveals 
a phase inconsistency of the TV* (1535) in the previous de- 
scription of dynamical generation of Ref. [Q . Part of the 
phase inconsistency could be traced back to the absence of 
the A^*(1650) resonance which strongly affects the proper- 
ties of the A^*(1535) through resonance interference. Thus, 
in an extension of the original model of Ref. [5] , we allow 
for two genuine 3-quark resonances: one to account for 
the A''*(1650), and another one to replace the dynamically 
generated Af*(1535), if the fit prefers this solution. 

In this work, we first show a fair agreement of the 
present model with the data close to threshold for n^p 
photoproduction. In the fit of the A^* region, one of the 
genuine poles indeed accounts for the A^*(1650). However, 
the other one, instead of replacing the dynamically gen- 
erated A^*(1535), moves far into the complex plane and 
provides an almost energy independent background; it ac- 
counts for background processes not explicitly included in 
the present photoproduction model such as t-channel vec- 
tor meson exchanges with anomalous photon couplings. 

In any case, the dynamically generated A'* (1535) pole 
appears as a stable object, resistant to changes of the fit. 
While its position changes significantly compared to the 
original model of Ref. [S] , due to an interference with the 
A* (1650), the present study shows that a dynamically 
generated Af*(1535), together with a genuine A^*(1650), 
can deliver a consistent picture simultaneously in nN — + 
nN and 7A^ nN. 

The dynamical generation is tied to the strong cou- 
plings to the KA and KS channels provided by the SU(3) 
Lagrangian. In this connection, we have also verified that 
the A* (1535) as a dynamical resonance disappears if the 
coupling strengths to the KA and KE channels are re- 
duced by about 40% - 50%. 

A detailed study of the analytic structure of the reac- 
tion amplitudes has revealed the role of the A^*(1535) and 
A* (1650) poles on other sheets; the pronounced cusp at 
the T]N threshold could be naturally explained. It appears 
as the result of an interplay of physical and hidden poles 
on different Riemann sheets. A virtual state in Sn below 
the nN threshold could be found that is quite stable and 
seems to be required by the sharp rise in Re that is 
seen in the partial wave analyses of, e.g., Ref. [JSj. It re- 
mains to be seen if this state is genuine or mocked up from 
t- and M-channel cuts that are not explicitly included in 
the present model [55] . 

Furthermore, we have found that the A^*(1535) pole 
disappears on some sheets; this implies a model-independent 
caveat: the use of sub-threshold resonances in phenomeno- 
logical analyses is questionable. Thus, it is important to 
consider all relevant Riemann sheets and to pay special 
attention to which parts of the physical axis they are con- 
nected; sub-threshold and above-threshold poles, that are 
easily overlooked, have important consequences for the 
physical amplitude and should be considered. 
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While the present model describes the S partial waves 
of the studied reactions well, there are residual discrepan- 
cies in the simultaneous description of the low and high 
energy regions. This could be traced back to the fact that 
the subtraction parameters are constants in energy; for a 
satisfactory description of the 5- wave amplitudes covering 
the energy region from threshold to second resonance re- 
gion, higher orders in the chiral meson-baryon interaction 
should be considered. Furthermore, although small for the 
5*11 amplitude, the ttttN channel should be included in the 
model. 

We have shown that the existing ttN data do not rule 
out the description of the A^*(1535) as a dynamically gen- 
erated resonance. However, whether or not this scenario 
is indeed the case, still remains to be seen. In particular, 
the influence of the strong couplings to the KA and KU 
channels -which are responsible for the dynamical genera- 
tion of this resonance- on the higher partial waves should 
be investigated as mentioned in the Introduction. 

While we have obtained a fair data description of the 
ttN final state in the present investigation with a dynami- 
cally generated iV*(1535), the next logical step is to study 
the other final states {rjN, KA, and KE) - already in- 
cluded in the model as intermediate states - in pion- and 
photon-induced reactions as well as in electroproduction 
to put further constraints on the model. 
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